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RESUMEN 

El resumen sera traducido al espanol por los editores. 
A new model for the postcollapse equilibrium structure of virialized objects which 
condense out of the cosmological background universe is described and compared 
with observations and simulations of cosmological halos. The model is based upon 
the assumption that virialized halos are isothermal, which leads to a prediction 
of a unique nonsingular isothermal sphere for the equilibrium structure, with a 
core density which is proportional to the mean background density at the epoch 
of collapse. These predicted nonsingular isothermal spheres are in good agreement 
with observations of the internal structure of dark-matter-dominated halos from 
dwarf galaxies to X-ray clusters. Our model also reproduces many of the average 
properties of halos in CDM simulations to good accuracy, suggesting that it is a 
useful analytical approximation for halos which form from realistic initial conditions. 
While N-body simulations find profiles with a central cusp, our nonsingular model 
matches the simulated halos outside the innermost region well. This model may 
also be of interest as a description of halos in nonstandard CDM models like self- 
interacting dark matter, which have been proposed to eliminate the discrepancy 
between the cuspy halos of standard CDM simulations and observed halos with 
uniform-density cores. 

ABSTRACT 

A new model for the postcollapse equilibrium structure of virialized objects 
which condense out of the cosmological background universe is described and com- 
pared with observations and simulations of cosmological halos. The model is based 
upon the assumption that virialized halos are isothermal, which leads to a predic- 
tion of a unique nonsingular isothermal sphere for the equilibrium structure, with 
a core density which is proportional to the mean background density at the epoch 
of collapse. These predicted nonsingular isothermal spheres are in good agreement 
with observations of the internal structure of dark-matter-dominated halos from 
dwarf galaxies to X-ray clusters. Our model also reproduces many of the average 
properties of halos in CDM simulations to good accuracy, suggesting that it is a 
useful analytical approximation for halos which form from realistic initial condi- 
tions. While N-body simulations find profiles with a central cusp, our nonsingular 
model matches the simulated halos outside the innermost region well. This model 
may also be of interest as a description of halos in nonstandard CDM models like 
self-interacting dark matter, which have been proposed to eliminate the discrep- 
ancy between the cuspy halos of standard CDM simulations and observed halos 
with uniform-density cores. 
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1. INTRODUCTION 

The question of what equilibrium structure results when a density perturbation collapses out of the expand- 
ing background universe and virializes is central to the theory of galaxy formation. The nonlinear outcome 
of the growth of Gaussian-random-noise cosmological density fluctuations due to gravitational instability in a 
hierarchical clustering model like CDM is not amenable to direct analytical solution, however. Instead, nu- 
merical simulations are required. As a guide to understanding these simulations, as a check on their accuracy, 
and as a means of extrapolating from simulation results of limited dynamic range, analytical approximations 
are nevertheless an essential tool. One such tool of great utility has been the solution of the spherical top-hat 
perturbation problem (cf. Gunn & Gott 1972, Padmanabhan 1993). As used in the Press-Schechter ("PS") 
approximation (Press & Schechter 1974) and its various refinements, the top-hat model serves to predict well 
the number density of virialized halos of different mass which form at different epochs in N-body simulations. 
An analytical model for the internal structure (e.g. mass profile, temperature, velocity dispersion, radius) of 
these virialized halos would be a further tool of great value for the semi-analytical modelling of galaxy and 
cluster formation, therefore. Here we shall summarize our attempt along these lines. 

Earlier work adopted crude approximations which used the virial theorem to match a collapsing top-hat 
perturbation either to a uniform sphere or to a singular isothermal sphere, with the same total energy as 
the top-hat. Our first motivation, therefore, is simply to improve upon this earlier treatment by finding a 
more realistic outcome for the top- hat problem. As a starting point, we shall adopt the assumption that 
the final equilibrium is spherical, isotropic, and isothermal, a reasonable first approximation to the N-body 
and gasdynamic simulation results of the CDM model. As we shall see, the postcollapse analytical solution 
derived from this assumption quantitatively reproduces many of the detailed properties of the halos found in 
those simulations, so we are encouraged to believe that our approximation is well justified. Our model is in 
disagreement, however, with the N-body simulation result that, in their very centers, dark-matter-dominated 
halos have cuspy profiles (e.g. Navarro, Frenk, & White 1997; "NFW"). By contrast, our model predicts a small, 
but uniform density core, as required to explain the observed dwarf galaxy rotation curves and cluster mass 
profiles inferred from gravitational lensing. This discrepancy between the cuspy profiles of the N-body results 
and the observed dark-matter-dominated halos has led recently to a reexamination of the cold, collisionless 
nature of CDM, itself, and the suggestion that a variation of the microphysical properties of the dark matter 
might make it more "collisional" , enabling it to relax dynamically inside these halos so as to eliminate the 
central cusp (e.g. Spcrgcl & Stcinhardt 1999). While the details of this suggestion are still uncertain, our 
model serves to predict its consequences, to the extent that we are able to ignore the details of the relaxation 
process inside the halo and approximate the final equilibrium as isothermal. In what follows, we shall describe 
our model and compare its predictions both with CDM simulation results and with observations of dwarf galaxy 
rotation curves and galaxy clusters. 

2. THE TRUNCATED ISOTHERMAL SPHERE MODEL 

Our model, as described in Shapiro, Iliev, and Raga (1999) for an Einstein-de Sitter universe and generalized 
to a low-density universe, either matter-dominated or a flat one with a positive cosmological constant, in Iliev 
& Shapiro (2000), is as follows: An initial top-hat density perturbation collapses and virializes, which leads 
to a truncated nonsingular isothermal sphere in hydrostatic equilibrium (TIS) , a solution of the Lane-Emden 
equation (appropriately modified in the A ^ case). Although the mass and total energy of the top- hat 
are conserved thru collapse and virialization, and the postcollapse temperature is set by the virial theorem 
(including the effect of a finite boundary pressure), the solution is not uniquely determined by these requirements 
alone. In order to find a unique solution, some additional information is required. We adopt the anzatz that 
the solution selected by nature will be the "minimum-energy solution" such that the boundary pressure is that 
for which the conserved top-hat energy is the minimum possible for an isothermal sphere of fixed mass within 
a finite truncation radius. As a check, we appeal to the details of the exact, self-similar, spherical, cosmological 
infall solution of Bcrtschinger (1985). In this solution, an initial overdensity causes a continuous sequence of 
spherical shells of cold matter, both pressure-free dark matter and baryonic fluid, centered on the overdensity, 
to slow their expansion, turn-around and recollapse. The baryonic infall is halted by a strong accretion shock 
while density caustics form in the collisionless dark matter, instead, due to shell-crossing. The postcollapse 



EQUILIBRIUM STRUCTURE OF COSMOLOGICAL HALOS 



3 








■1.5 — 1 — 0.5 
l°gio(i'/ r „) 



/ 2 GMm \ 
\5 r vir ) 



kg T v j r 
' 1 GM 

Po/Pt 

(P)/Pt 

r t /r 

(P) 

Pcrit (^coll) 



sus 

0.5 



Tabic 1 
SIS 
0.417 

3 



TIS (CI = 1) 
0.554 

2.16 



1 


oo 


514 


1 


3 


3.73 


NA - 


oo 


29.4 


18tt 2 


"-(!)' 


130.5 



Fig. 1. Density profile of truncated isothermal sphere which forms from the virialization of a top-hat density 
perturbation in a matter-dominated universe. Radius r is in units of r m - the top-hat radius at maximum 
expansion, while density p is in terms of the density psus of the standard uniform sphere approximation for 
the virialized, post-collapse top-hat. Bottom panel shows logarithmic slope of density profile. 



virializcd object we wish to model is then identified with that particular shock- and caustic-bounded sphere 
in this infall solution for which the mass and total energy are the same as those of our top-hat before collapse 
and the trajectory of its outermost mass shell was identical to that of the outer boundary of our collapsing 
top-hat at all times until it encountered the shock. This spherical region of post-shock gas and shell-crossing 
dark matter in the infall solution is very close to hydrostatic and isothermal and has virtually the same radius 
as that of the minimum-energy solution for the TIS. This confirms our "minimum-energy" anzatz and explains 
the dynamical origin of the boundary pressure implied by that solution as that which results from thcrmalizing 
the kinetic energy of infall. 

With this "minimum-energy" anzatz, we find that a top- hat perturbation collapse leads to a unique, non- 
singular TIS, which yields a universal, self-similar density profile for the postcollapsc equilibrium of cosmic 
structure. Our solution has a unique length scale and amplitude set by the top-hat mass and collapse epoch, 
with a density proportional to the background density at that epoch. The density profiles for gas and dark 
matter are assumed to be the same (no bias). The final virialized halo has a flat density core. 

Case I: matter-dominated cases, both flat and low density (see Fig. 1). The core size is r$ = 0.034 x 
radius r t , where rt is the size of the halo (i.e. truncation radius). The central density is po = 514 x surface 
density p t . The ID velocity dispersion oy of the dark matter and the gas temperature T are then given by 
Cy — ksT/ (prrip) — 47rGpo r o- [Note: this ro = ^o, King/3, where ro,King is the core radius defined as the "King 
radius" by Binney & Trcmaine (1987, equ. [4-124b])]. Compared to the standard uniform sphere (SUS) and 
singular isothermal sphere (SIS) approximations (see Table 1), the temperature is T = 2.16 Tsus = 0.72 Tsis- 
At intermediate radii, p drops faster than r~ 2 . 

Case II: flat, A ^ models. The profile varies with epoch of collapse, approaching the universal shape 
of case I above for early collapse. For example, for Cl = 1 — A = 0.3 and z co \\ = (0;0.5;1), we obtain 
r t /r = (30.04; 29.68; 29.54), p /p t = (529.9; 520.8; 517.2), and T/T sus = (2.188; 2.170; 2.163), respectively. 



3. DWARF GALAXY ROTATION CURVES AND THE V M ax - Umax CORRELATIONS 

The TIS profile matches the observed mass profiles of dark-matter-dominated dwarf galaxies. 

The observed rotation curves of dwarf galaxies can be fit according to the following density profile with a finite 
density core (Burkert 1995): 

/ \ PO, Burkert / -. \ 

P[T) ~ (r/r c + l)(r 2 /r 2 + l)' 1 ' 

The TIS profile gives a nearly perfect fit to the Burkert profile, with best-fit parameters po. Burkert/ Po.tis — 1-216, 
r c /r ,Tis = 3.134 (see Fig. 2a). 
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Fig. 2. (a) (left) Rotation Curve Fit. Solid line = Best-fit TIS; Dashed line = Burkcrt profile, where <j\ is = 
(v 2 }/3 = feT/m. (b) (right) fmax-rmax correlations. Observed dwarf galaxies (triangles) and LSB galaxies 
(squares) from Kravtsov et al. (1998); Burkert: fit to data (Burkcrt 1995); SCDM: fi = 1, A = 0, a 8h -i = 0.5 
(cluster-normalized); OCDM: Hq = 0.3, Ao = (all COBE-normalized, subscript indicates n-value of tilt); 
ACDM: Q = 0.3, A = 0.7, (COBE-normalized; no tilt); h = 0.65 for all. 



How well does this best- fit TIS profile predict v maxi the maximum rotation velocity, and the radius , r max , 
at which it occurs in the Burkert profile? We find r max burkert / r max,T is — 1-13, v max . Burkert / v max,T is — 

1.01 

(i.e. excellent agreement). 

The TIS halo model explains the observed correlation of v max and r max for dwarf spiral and 
LSB galaxies, when the TIS halo model is combined with the Press-Schcchter model to predict the typical 
collapse epoch for objects of a given mass (i.e. the mass of the 1-a fluctuations vs. z co u) (See Fig. 2b). Both 
of the flat, untilted CDM models plotted, cluster-normalized Einstcin-de Sitter and COBE-normalized, flat, 
low-density models (Oo = 0.3 and Ao = 0.7), as well as the slightly tilted (n = 1.14) open model (f2o = 0.3) 
yield a reasonable agreement with the observed v max — r max relation, while the untilted (n = 1) and strongly 
tilted (n — 1.3) open models do not agree with the data. 

4. GALAXY HALO M - a v RELATION 

Our TIS halo model predicts the velocity dispersion of galactic halos of different mass which 
form in the CDM model according to N-body simulations. Antonuccio-Delogu, Becciani, & Pagliaro 
(1999) used an N-body treecode at high-res (256 3 particles) to simulate galactic halos in regions of a single 
and of a double cluster. They found that the agreement with the TIS model is quite good, much better than 
with either of the other two models they considered, namely the singular isothermal sphere and the peak-patch 
model of Bond & Myers (1996). 

5. COMPARISONS WITH GALAXY CLUSTER OBSERVATIONS AND SIMULATIONS 

The TIS halo model predicts the internal structure of X-ray clusters found by gas-dynamical/N- 
body simulations of cluster formation in the CDM model. Our TIS model predictions agree astonish- 
ingly well with the mass-temperature and the radius-temperature virial relations and integrated mass profiles 
derived from numerical simulations by Evrard, Metzler and Navarro (1996; "EMN"). Apparently, these simula- 
tion results are not sensitive to the discrepancy between our prediction of a finite density core and the N-body 
predictions of a density cusp for clusters in CDM. Let X be the average overdensity inside the sphere of radius 

1 /2 

r, X = (p(r))/pi,. Then the radius-temperature virial relation is defined as rx = ri (A)(T/10keV) ' Mpc, 
and the mass-temperature virial relation by Mx = Mi (A)(T/10keV) 1,/2 /i _1 10 15 M Q . A comparison between 
our predictions of the mass - temperature relation r w (X) and the results of EMN is given in Fig. 3a. For 
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Fig. 3. (a) (left) Cluster radius-temperature virial relation (z = 0). (triangles) CDM simulation results as fit 
by Evrard, et al. (1996); (continuous line) TIS prediction, (b) (right) Projected surface density of cluster CL 
0024+1654 inferred from lensing measurements, together with the best-fit TIS model. 



the mass-temperature virial relation EMN obtain Mi (500) = 1.11 ± 0.16 and Mi (200) = 1.45, while our TIS 
solution yields Mi (500) = 1.11 and Mi (200) = 1.55, respectively. 

The TIS model for the internal structure of X-ray clusters predicts gas density profiles p gaa (r) 
and X-ray brightness profiles 1(9) which are well-fit by the standard (3 profile, 



Pgas 



Po 



(l + r 2 /rl) 



3/3/2 • 



(i + e 2 /e, 



2 //)2\3/3-l/2^ 



(2) 



with /3-values for the TIS /3-fit which are quite close to those of simulated clusters in the CDM model but 
somewhat larger than the conventional observational result that (3 ~ 2/3 (see tables below). However, recent 
X-ray results suggest that the true (3- values are larger than 2/3 when measurements at larger radii are used 
and when central cooling flows are excluded from the fit. 
I(r) (observations) (3 



Jones and Foreman (1999) 
Jones and Foreman (1992) 
Balland and Blanchard (1997) 

Durret et al. (2000) 

Vikhlinin, et al. (1999) 
(fit by Henry 2000) 



0.4-0.8, ave. 0.6 

- 2/3 
0.57 (Perseus) 
0.75 (Coma) 
0.53 (incl. cooling flow) 
0.82 (excl. cooling flow) 
0.7-0.8 



Pgas(^) (simulations) 



Metzler and Evrard (1997) 

Eke, Navarro, and Frenk (1998) 
Lewis et al. (1999) (adiabatic) 
Takizawa and Mincshige (1998) 
Navarro, Frenk, and White (1995) 



0.826 (DM) 
0.870 (gas) 
0.82 
~ 1 
- 0.9 
0.8 



TIS /3-fit (r c /r , T /s = 2.639) 



0.904 



TIS /J-fit {r c /r ,Tis = 2.416) 



0.846 



The TIS halo model can explain the mass profile with a fiat density core measured by Tyson, 
Kochanski, and DelPAntonio (1998) for cluster CL 0024+1654 at z = 0.39, using the strong 
gravitational lensing of background galaxies by the cluster to infer the cluster mass distribution. 

The TIS model not only provides a good fit to the shape of the projected surface mass density distribution of 
this cluster within the arcs (see Fig. 3b), but when we match the central value as well as the shape, our model 
predicts the overall mass, and a cluster velocity dispersion in close agreement with the value a v — 1200 km/s 
measured by Dressier and Gunn (1992). By contrast, the NFW fit which Broadhurst et al. (2000) reports can 
model the lensing data without a uniform-density core predicts oy much larger than observed. 
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6. SUMMARY 



• The TIS profile fits dwarf galaxy rotation curves; combined with the Press-Schechter approximation, it 
predicts the observed v ma x — r max relation for dwarf and LSB galaxies. 

• The TIS predicted M — ay relation agrees with high resolution N-body simulations of galactic halo 
formation by Antonuccio-Delogu et al. (1999). 

• The predicted mass-radius-temperature scaling relations and integrated mass profile of the TIS model 
match simulation results for clusters in the CDM model in detail. Our solution derives the empirical 
fitting formulae of Evrard, et al. (1996), which also agree well with X-ray cluster observations at z = 0. 

• The TIS X-ray brightness profile matches the /3-fit profile with (3 m 0.9, larger than typically reported by 
X-ray observers, but very close to the results of gas-dynamical/N-body simulations of X-ray clusters in 
the CDM model. 

• The TIS solution fits the cluster mass profile with uniform-density core derived from strong gravitational 
lensing data by Tyson et al. (1998) for CL 0024+1654 within the arcs, while accurately predicting the 
observed av on larger scales, too. 

This work was supported by grants NASA NAG5-2785, NAG5-7363, and NAG5-7821, NSF ASC-9504046, 
and Texas ARP 3658-0624-1999, and a 1997 CONACyT National Chair of Excellence at UNAM for PRS. 
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